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ABSTRACT 

We identify a large sample of isolated bright galaxies and their fainter satellites in 
the 2dF Galaxy Redshift Survey (2dFGRS). We analyse the dynamics of ensembles 
of these galaxies selected according to luminosity and morphological type by stacking 
the positions of their satellites and estimating the velocity dispersion of the com- 
bined set. We test our methodology using realistic mock catalogues constructed from 
cosmological simulations. The method returns an unbiased estimate of the velocity 
dispersion provided that the isolation criterion is strict enough to avoid contami- 
nation and that the scatter in halo mass at fixed primary luminosity is small. Us- 
ing a maximum likelihood estimator that accounts for interlopers, we determine the 
satellite velocity dispersion within a projected radius of 175 h^^ kpc. The disper- 
sion increases with the luminosity of the primary and is larger for elliptical galax- 
ies than for spiral galaxies of similar bj luminosity. Calibrating the mass-velocity 
dispersion relation using our mock catalogues, we find a dynamical mass within 
175 /i-i kpc of Mi75/h-iM0 ~ x {l^^Jl^*) for elliptical galaxies and 

Mi75/h-iM0 ~ 6.3t|5 X 10" (Lb,j/L*)i-6 for spiral galaxies. Finally, we compare our 
results with recent studies and investigate their limitations using our mock catalogues. 

Key words: galaxies: halos - galaxies: kinematics & dynamics - galaxies: mass - 
galaxies: satellites - galaxies: spiral - surveys 



,;_( 1 INTRODUCTION 



The view that galaxies are surrounded by large dark mat- 
ter halos dates back more than 30 years to the pioneering 
study of the rotation curve of M32 by Rubin & Ford (1970). 
Extended galactic halos are, in fact, a generic feature of the 
cold dark matter model of galaxy formation (Blumenthal et 
al 1984, Frenk et al 1985), but this fundamental theoretical 
prediction has limited observational support. Zaritsky et al. 
(1993) attempted to measure the mass and extent of dark 
matter halos by analysing the dynamics of satellite galax- 
ies found around 'isolated' galaxies. Since galaxies generally 
have only a few detectable satellites, they used a method 
that consists of stacking satellites in a sample of primaries 
of similar luminosity. In spite of the small size of their rel- 
atively inhomogeneous sample, Zaritsky et al. (1993) were 
able to detect massive halos around isolated spiral galax- 
ies extending to many optical radii. Having nearly doubled 
their satellite sample to 115 members, Zaritsky et al. (1997b) 
confirmed their earlier claims including a puzzling lack of 
correlation between the velocity dispersion of the (stacked) 
satellite system and the luminosity of the primary. 
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More recently, McKay et al. (2002) performed a similar 
analysis on data from the Sloan Digital Sky Survey (SDSS; 
York et al. 2000). They compared mass estimates derived 
from satellite dynamics to those derived from weak lens- 
ing analyses of the same data (McKay et al. 2001). With 
a much larger sample than that of Zaritsky et al (1997), 
they were able to detect a correlation between satellite ve- 
locity dispersion and primary luminosity. This trend was 
confirmed by Prada et al. (2003) who also used SDSS data. 
Although they are both based on SDSS data, these two stud- 
ies find results that, while consistent at first sight, are in, 
fact, somewhat contradictory. For example, although Prada 
et al (2003) found a strong dependence of satellite velocity 
dispersion on galactrocentric distance, their measured veloc- 
ity dispersion within a radius of 125 kpc is similar to 
the values obtained by McKay et al. (2002) at a radius of 
275 kpc. Discrepant results were also found by Brain- 
erd & Specian (2003) who applied the same technique to 
the early, "100k" data release of the 2degree-Field Galaxy 
Redshift Survey (2dFGRS; CoUess 2001) and derived sateb 
lite velocity dispersions which are in qualitative and quan- 
titative disagreement with those of Zaritsky et al. (1997), 
McKay et al. (2002) and Prada et al. (2003). A more 
extensive analysis of the complete 2dFGRS (CoUess et al. 
2003) by Brainerd (2005) also led to disagreements with the 
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results of earlier work. This somewhat confused picture of 
satellite dynamics is due in large part to different choices of 
primary galaxy samples and to differences in the modelling 
and analysis methods. 

This paper has multiple aims. Firstly, we carry out a 
new analysis of the dynamics of satellites around bright 
galaxies of different morphological types selected from the 
full 2dFGRS. The goal is to constrain the velocity dispersion 
and mass of their dark matter halos and we therefore select 
a sample of isolated galaxies chosen according to strict cri- 
teria. Secondly, we investigate the reliability and accuracy 
of commonly used dynamical analysis methods. For this, we 
make extensive use of realistic mock catalogues constructed 
from large cosmological N-body simulations and different 
semi-analytic galaxy formation models (Cole et al. 2000; 
Springel et al. 2001). A similar approach, but in a differ- 
ent context, was adopted by van den Bosch et al. (2004). 
Finally, we attempt to understand the root cause of the dif- 
ferences found in previous work, again relying on the use of 
realistic mock catalogues. 

The paper is organised as follows. In Section (2] we 
briefly present some of the characteristics of the 2dFGRS 
data and simulations used for our analysis. In Section |31 
we describe the satellite sample selection scheme, together 
with its robustness to changes in the selection parameters. 
The analysis of the 'stacked' satellite velocity distribution 
is carried in Section |4] while, in Section [5] we present ve- 
locity dispersion estimates for our mock catalogues and for 
2dFGRS primaries split according to luminosity and mor- 
phological type. Using a model for the relationship between 
dark halo mass and satellite velocity dispersion, we give, 
in Section [S] an estimate of the mass of the halos around 
2dFGRS galaxies. In Section [T] we compare our results with 
those of previous studies, and we conclude in Section [S] 

2 THE DATA 

2.1 The 2dFGRS data 

Detailed descriptions of the construction of the 2dF Galaxy 
Redshift Survey (2dFGRS) and its properties are given in 
CoUess et al. (2001; 2003). In summary, galaxies are se- 
lected down to a magnitude limit of bj ~ 19.45 from the 
full 2dFGRS catalogue, with ~ 225 000 galaxies having a 
"good quality" redshift measurement. We restrict our anal- 
ysis to the two large contiguous volumes of the survey, one 
centred on the Southern Galactic Pole (SGP) and the other 
close to the direction of the Northern Galactic Pole (NGP). 

Three limitations of the 2dFGRS catalogue are relevant 
for this work. Firstly, the 2dFGRS source catalogue, based 
on the APM galaxy catalogue, is not complete. By compar- 
ing the 2dFGRS with the SDSS, Norberg et al. (2002) esti- 
mated the completeness of the 2dFGRS to be ~ 91% ± 2%, 
and ascribed the incompleteness primarily to misclassifica- 
tion of APM images. Misclassification of close galaxy pairs 
will cause some true pairs to be missed from our sample 
of primaries. The 2dFGRS suffers from an additional form 
of close pair incompleteness due to 'fibre collisions' during 
the spectroscopic observations. The observing strategy em- 
ployed in the 2dFGRS, consisting of a set of overlapping tiles 
which are successively observed, was designed to minimise 
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Figure 1. A schematic of the isolation criterion. The figure on the 
left shows the projection on the sky of the cylinder that defines 
the primary isolation criterion. The figure on the right shows the 
configuration as seen along the redshift axis. The central point in 
each panel corresponds to a primary; only galaxies located in the 
shaded area are considered as potential satellites. 

the number of 'fibre collisions', and the remaining incom- 
pleteness is very precisely quantified (CoUess et al. 2001). 
Finally, the rms accuracy of redshift measurements for a 
typical galaxy is 85 km s~^ and tends to be slightly larger 
for the faintest galaxies and slightly smaller for the brighter 
ones (CoUess et al. 2001). 

Since we are primarily interested in the velocity dis- 
persion at large radius, the close pair incompleteness in the 
catalogue does not have an important effect on our conclu- 
sions. Nevertheless, it is important to model carefully both 
the incompleteness and the velocity errors and to include 
these in the construction of our mock galaxy catalogues. 

Finally, we make use of the 'eyeball' morphological clas- 
sification carried out by Loveday et al. (1996) based on the 
APM images. This is available for 80% of the central galax- 
ies in our sample. We prefer this classification to the more 
objective spectral classification of Madgwick et al. (2002) 
because aperture effects are important for our sample of rel- 
atively nearby galaxies. However, we have repeated the anal- 
ysis of section using subsamples defined by their spectral 
classification and find no difference in our results within the 
errors. 



2.2 The ACDM Simulation 

We use mock 2dFGRS catalogues constructed from cosmo- 
logical simulations in order to assess the extent to which the 
various limitations of the data affect our results. In particu- 
lar, we use the mocks to investigate possible systematic ef- 
fects arising from our method for selecting satellites, as well 
as from our procedure of stacking satellites together. The 
mock catalogues allow us also to investigate the effects of 
redshift space distortions and redshift measurement errors. 
In the simulations we, of course, know not only the redshifts 
of galaxies but also their distances. In what follows, we use 
the term "real space" to refer to measurements that make 
use of the true 3-d position and the term "redshift space" 
to refer to measurements that make use of pseudo 3-d posi- 
tions, i.e. those for which the distance to the galaxy is given 
by the sum of the pure Hubble fiow distance and the peculiar 
velocity along the line-of-sight, in units of /i^^Mpc. 

To construct mock 2dFGRS catalogues, we use a high- 
resolution N-body simulation of a flat, A-dominated cold 
dark matter universe with the following parameters: matter 
density, Qrn ~ 0.3; cosmological constant term, Qa = 0.7; 
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Table 1. Properties of the combined NGP & SGP satellite samples around bright galaxies, for different values of N^i^l a-nd for different 
primary morphological type. The numbers quoted in brackets are after small groups are excluded from the full satellite sample, rp is 
in kpc. 



Sample Primary Type N^^^i A^prim A'sat A^sat (fint < >" < »'cxt ) 

< rp < 175 40 < rp < 175 



2dFGRS 


any 





362 


(357) 


642 


(588) 


273 


241 


2dFGRS 


any 


4 


571 


(564) 


1003 


(918) 


434 


383 


2dFGRS 


Spiral- Irregular 





203 


(203) 


322 


(322) 


141 


120 


2dFGRS 


Spiral-Irregular 


4 


311 


(311) 


465 


(465) 


210 


181 


2dFGRS 


EUiptical-SO 





85 


( 81) 


203 


(161) 


75 


69 


2dFGRS 


EUiptical-SO 


4 


141 


(135) 


338 


(265) 


118 


111 


Mock 


any 





387 


(387) 


648 


(648) 


384 


231 


Mock 


any 


4 


736 


(736) 


1226 


(1226) 


723 


442 



Hubble constant, Ho = 70 km s~^Mpc~^; index of primor- 
dial fluctuation power spectrum, n = 1; and present-day 
fluctuation amplitude as — 0.9. The simulation followed 
400'^ particles in a box of of side 110 /i~^Mpc (see Gao et al. 
(2004) for a full description of the simulation). The simula- 
tion was populated with galaxies by applying the "Munich" 
semi-analytic model of galaxy formation to the merger trees 
of each halo (Springel et al. 2001). In this model, galaxies re- 
side in resolved halos and their subhalos. When a subhalo is 
no longer resolved, the galaxy is placed on the most bound 
particle of the subhalo when it was last resolved and an 
analytic dynamical friction calculation is used to determine 
when a satellite merges with the central galaxy. The free pa- 
rameters of the model are tuned to match the TuUy-Fisher 
relation, the B-band cluster galaxy luminosity function and 
the overall two-point correlation function (Springel et al. 
2001). An interesting feature of this semi-analytic galaxy 
formation model is the generation of a velocity bias between 
galaxies and dark matter, as function of halo radius. 

In order to obtain as close a match as possible between 
the 2dFGRS selection function and that in the mocks, we 
rescale the luminosities of the model galaxies preserving the 
luminosity ranking so that their luminosity function exactly 
matches that of the 2dFGRS (Norberg et al. 2002). The 
required rescaling can be as large as a magnitude for some 
of the mock galaxies, but the differential rescaling for galax- 
ies brighter than Mbj — 51ogioh< —18, which make up the 
bulk of our sample of primaries, is small. Thus, for these 
brighter galaxies, the magnitude differences are quite accu- 
rately preserved. We then extract magnitude-limited cata- 
logues of galaxies to the same magnitude limit and with the 
same geometry as the real 2dFGRS, as described by Norberg 
et al. (2002). This requires using three periodic replications 
of the simulation cube. (We have checked that removing any 
duplicated systems does not affect any or our results.) 

The last step is to create 'sampled' mocks by applying 
the 2dFGRS masks (including the redshift incompleteness 
mask) and the 2dFGRS photometric errors, as explained in 
greater detail in Norberg et al. (2002). Finally, we add to 
each observed velocity an "observational error" randomly 
sampled from a Gaussian of width 80 km s~^. This value is 
a compromise between the Gaussian errors measured from 
repeat observations for our sample of 2dFGRS primaries 
(o- ~ 70 ± 7 km s"^) and satellites (cr ~ 86 ± 5 km s"^). 



3 SATELLITE SAMPLE 

We begin this section by explaining the method used to de- 
fine the satellite sample, which is applied to both the 2dF- 
GRS data and the mock catalogues. We then briefly consider 
the robustness of the satellite properties to variations of the 
selection parameters, an issue that we address further in the 
appendices. Finally, we present some general properties of 
the satellite samples used in this paper. 

3.1 Satellites around Isolated Primaries 

Since the main purpose of our analysis is to constrain the 
mass of the galactic halo using the dynamics of satellites, we 
require a sample of isolated primaries. To construct it, we 
begin by excluding regions in the 2dFGRS that could be con- 
taminated by large clusters. Specifically, we exclude regions 
lying within three projected Abell radii (3 x 1.5^^^ kpc) 
and 3000 km s~^ of the centre of clusters in the catalogue of 
Dalton et al. (1997). Next, we select a sample of bright, iso- 
lated primaries by requiring that they satisfy the following 
criteria: 

• the local 2dFGRS magnitude limit should be at least 
SrUs = 2.2 fainter than the primary 

• all the neighbouring galaxies within AV — IV^p"^™ — 
ygai| ^ 2400 km s~^ and within a projected radius Srs < 
400 kpc should be faint enough to satisfy foj''' — bjj"^'™ > 
StUs = 2.2 

• all neighbouring galaxies within AV — jV''"^'™ — 
ygai| ^ 2400 km s~^ and within a projected radius Srt < 
1000 kpc should satisfy 6f ' - 6f ™ > Smi, = 0.8, 

where the projected radius between two galaxies, at posi- 
tions ri and r^, is defined by: 

e, „iri + r2i /l -cos(a) . ri • r2 , 

da — 2 \ / with cos(a) — - — - — - . 

2 V H-cos(a) ^ ' |ri| |r2| ^ 

All galaxies that lie within a projected distance < 
400 h^^ kpc of a primary and have relative velocity AVa = 
I y prim _ y^gai| ^ ^jOO km s'^ are considered as potential 
satellites. The isolation criterion is presented in schematic 
form in Fig. [T] 

Our adopted value of 5ms corresponds to a factor of 
8 in luminosity and a similar factor in mass. The primary 
motivation behind this choice is to ensure that the satellites 
are small enough to produce only minor perturbations in the 
gravitational potential of the system. As shown in section O 
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our galaxy mock catalogues indicate that our adopted value 
is adequate. A more detailed discussion of the precise choice 
of StUs may be found in Appendix A. 

Since not all galaxies in the 2dFGRS have a measured 
redshift, the primary isolation criterion could be violated by 
galaxies that lack a redshift measurement. We can guard 
against this by eliminating all primaries that could have 
their isolation criterion violated by such galaxies. A less 
restrictive condition is to accept only those for which it 
could be violated by at most Nvioi galaxies. A conservative 
estimate of Nvioi follows from taking all galaxies without 
redshift to be at the redshift of the primary and checking 
whether this would cause the primary to violate the isola- 
tion criterion. In what follows, we adopt a value of Nvioi = 4. 
We have checked that none of our results are influenced by 
the precise choice of Nvioi used. The only effect of adopting 
Nvioi = 4, rather than Nvioi = 0, is to increase the number of 
satellites and hence the signal-to-noise of our measurements. 
With this choice of Nvioi, the satellite sample is ~ 55 % larger 
than for Nvioi = (~ 45 % for spirals and ~ 65 % larger 
for ellipticals). The largest increase occurs for the faintest 
absolute magnitude bins for which the number of satellites 
nearly doubles. 

With this algorithm and using the values of the selec- 
tion parameters specified above, we identify 571 primary 
galaxies surrounded by 1003 satellites. The algorithm also 
detects over 1500 isolated galaxies without spectroscopically 
confirmed satellites brighter than the local magnitude limit, 
but these have over 2200 neighbours without measured red- 
shift which could in principle be satellites. Some of these 
statistics are summarised in Table [1] 

Applying the same selection criterion to the mock cat- 
alogues, we identify, in real space, 750 primary galaxies sur- 
rounded by 1241 satellites. In real space, the isolation crite- 
rion is slightly different: instead of two cylinders of length 
AV and AVs, we use two spheres of radius Sri, and Svs re- 
spectively. The requirements on the magnitude differences 
remain the same. All galaxies within Srs of the primary 
are considered as satellites. In redshift space, including ve- 
locity errors, we find 736 primary galaxies surrounded by 
1226 satellites. As for the real data, the algorithm detects a 
further 1500 primaries without any confirmed satellite, but 
with 1700 possible candidates without measured redshifts. 

3.2 Robustness of the satellite detection 
algorithm 

It is important to test the dependency of the satellite de- 
tection algorithm on the values of the selection parameters. 
We find that the size of the satellite sample and its proper- 
ties are not very sensitive to the specific values of the inner 
and outer projected radii, so long as the exclusion crite- 
rion does not become too restrictive. The adopted values for 
these parameters represent a compromise between having a 
dynamically isolated system and a large sample of satellites. 

On the other hand, the inner and outer cylinder depths 
have a non-trivial influence on the satellite sample. First, if 
the 'velocity difference' between AVs and AV is less than 
then there is the potential risk of finding a single satel- 

1 i.e. \AV-AVs\ <AVs 
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Figure 2. Transverse distributions of satellites; comparison be- 
tween satellite samples extracted from the 2dFGRS (solid bold 
line), from the mocks in redshift space (dashed line) and in real 
space (dotted line). The vertical dashed line indicates the limiting 
projected radii, as used for the 'rp cut' (see text). The errorbars 
plotted assume Poisson statistics. 

lite galaxy associated with two difi'erent primaries. In such a 
case, neither primary can be considered isolated. Therefore 
in order to avoid this problem and obtain a self consistent 
isolation criterion, we always impose jAV — A14| > AVs. We 
note that none of the previous work in this subject has in- 
cluded such a constraint. It is not clear if their samples had 
any shared satellites and if they did how they were treated, 
but we note that for both the 2dFGRS data and our mock 
galaxy catalogues, this subtle problem does occur. 

Secondly, if the outer depth of the cylinder, i.e. AV, is 
very large, the isolation criterion is more stringent , and our 
catalogue will contain fewer primaries. Conversely if AV is 
too small the isolation criterion becomes too relaxed and we 
risk including non-isolated systems in our sample. 

Finally, if AVs increases the contamination of the satel- 
lite sample by interloper^ is increased. For instance, when 
increasing AVs from 600 km s~^ to 1800 km s~^, we see 
a flattening of the velocity distribution, which we interpret 
as being due to interlopers in our satellite catalogue. Con- 
versely too small a value of AVs would mean that the full 
width of the velocity distribution that we are trying to char- 
acterise would not be sampled. It is essential that AVs be 
greater than 3 to 4 ct, where a is rms width of the underlying 
satellite velocity distribution. 

Therefore the choice of AV = 2400 km s"^ and AVs = 

^ In the mock catalogues we can label any galaxy that is selected 
as a satellite, but does not reside within the dark matter halo 
of the primary galaxy as an interloper. In the real data this is 
not possible, but one can still statistically estimate the fraction 
of interlopers by their effect on the satellite velocity distribution 
(see section |4]2]l . 
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Figure 3. Statistics of primary galaxies and their associated satellites: in all panels we adopt a thick line style for 2dFGRS data and a thin 
line style for final mock data. All errorbars assume Poisson statistics. The top left panel shows the distribution of absolute magnitudes of 
satellites (solid line) and primaries (dot-dashed line) . The top right panel shows the distribution of the number of satellites per primary. 
The bottom left panel shows the distribution of the magnitude difference between primary and satellite. The bottom right panel shows 
the fraction of satellites per primary as function of rcdshift. All panels are done using the satellite samples obtained by applying the 
isolation criterion, rejecting small groups and with 'rp cut'. See text for more details. 



1200 km is a compromise between reducing the cata- 
logue contamination from interlopers, increasing the size of 
the satellite catalogue and allowing a robust velocity dis- 
persion estimate to be obtained. We note that the choice of 
the depth of the cylinders could be tuned with the size of 
the system considered to increase the efficiency with which 
'small' systems are detected. This is a point to which we 
return in section [S] 

In the appendices, we discuss in detail the isolation cri- 



teria proposed by McKay et al. (2002), Prada et al. (2003), 
Brainerd & Specian (2003) , van den Bosch et al. (2004) and 
Brainerd (2005). The summary of those findings is given in 
section [T] 



3.3 General Properties of Satellite Samples 

Before performing a detailed dynamical study, we focus 
briefly on some general properties of the satellite samples. 
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This leads us to make some additional cuts to improve the 
match between observed and mock samples. 

3. 3. 1 Transverse distributions of satellites 

In Fig. O we present the distribution of transverse separa- 
tions for three satellite samples. The mock satellite samples 
selected both in real and redshift space have centrally peaked 
distributions with the redshift space selected sample having 
the flatter distribution at separations rp > 100 kpc. In 
contrast the data has a almost fiat distribution throughout 
100 < Tp < 350 kpc, but with a significant drop in the 
central region. The reasons behind this difference on small 
scales are multiple: first of all, the 2dFGRS input catalogue 
lacks close galaxy pairs (Norberg et al. 2002) and this deficit 
is enhanced by the targets that are rejected due to fibre col- 
lisions; secondly galaxies are not point source objects, but 
extended objects on the sky, which means that the inner- 
most radial bin can suffer from projection effects which are 
not taken into account in the mocks (see van den Bosch et 
al. (2005) for a more detailed study of this particular is- 
sue). Finally, the large number of satellites in the innermost 
radial bin is not a generic model prediction: their number 
is sensitive to details of the dynamical friction prescription 
used in the semi-analytic model. 

We could model this limitation in the mock cata- 
logues, by e.g. implementing a supplementary incomplete- 
ness around each primary. However, we preferred not to 
add an arbitrary incompleteness model to our analysis, and 
hence choose to restrict some comparisons to satellite sam- 
ples with satellites rp > 40 kpc. Hereafter we refer to 
this as the 'rp cut'. The vertical line in Fig. [2] indicates this 
inner limiting projected radius. Beyond this radius, the dis- 
tributions of transverse separations for satellites from the 
mocks (in redshift space) and from the 2dFGRS are in ap- 
proximate agreement. 

3.3.2 Satellite sample properties 

In Fig. [3l we present general properties from the mock and 
data satellite samples (with light and bold line-styles respec- 
tively). The data and the semi-analytic mock satellite sam- 
ples have very similar properties: the peak of the satellite 
and primary absolute magnitude distributions are roughly 
the same for both samples, with primaries being typically 
one magnitude brighter than M* (top left panel of Fig. [S}; 
the fraction of satellites per primary does not vary signif- 
icantly as function of redshift; the overall distribution of 
satellite-primary magnitude differences is similar for both 
samples, as is the shape of the distribution of the number of 
satellites per primary (top right panel of Fig. |3]). 

The good match between data and mock satellite sam- 
ples, shown in Fig. [31 is achieved after applying both the 
'rp cut' and an additional cut to remove small groups from 
the 2dFGRS sample. We have removed all primaries which 
have 9 or more satellites. In the mock catalogue there is only 
one primary with more than 6 satellites, but in the 2dFGRS 
there are a few primaries, satisfying the isolation criterion, 
with 9 or more satellites. Most of them, more than 85%, 
have primaries that are elliptical galaxies. Removing these 
systems reduces our 2dFGRS satellite sample by 85 satellites 
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Figure 4. Velocity difference of satellite galaxies and primaries 
versus transverse separation for the 2dFGRS sample (bottom 
right panel) and for the mock samples: real space (top left panel); 
redshift space without velocity errors (bottom left panel) and with 
velocity errors (top right panel). The solid lines show respectively 
the 16"", 50"" and 84"^ percentiles of the distribution, in bins of 
40 kpc. This plot has not been corrected for interlopers. 

and reduces the the average number of satellites per primary 
from ~ 1.75 to ~ 1.63. The reason for removing them is to 
both achieve a better match between the mock and 2dF- 
GRS samples and to avoid our dynamical estimates being 
dominated by these small groups. We could have achieved 
this second goal by retaining the groups, but down weight- 
ing them by giving equal weight per primary rather than per 
satellite. However, we find that within the errors this does 
not change our results. 

In summary, with these extra restrictions, we end up 
with mock satellite samples which are rather similar to the 
2dFGRS satellite samples. 



4 MODELLING THE SATELLITE VELOCITY 
DISPERSION 

In this section, we consider the dynamical properties of the 
satellite samples obtained in section [3] We start by looking 
at the satellite velocity distributions, then address the issue 
of interlopers and background subtraction. Once they are 
well understood, we devise a method to estimate the satellite 
velocity dispersion of stacked primaries. 

4.1 Velocity Distribution of Satellites 

Fig. 13] is a scatter plot showing, for the 2dFGRS data and 
the mock samples, the satellite galaxy velocity difference 
(with respect to its associated primary) versus its projected 
distance from the primary. In all panels, the velocity dis- 
tribution is rather symmetric around zero for all projected 
distances. The precise choice of the cylinder depth (as fixed 
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Figure 5. Velocity distributions of satellites within a projected 
radius of 400 kpc: comparison between 2dFGRS satellite 
sample (solid bold line) and the mock samples in real and redshift 
space (with and without 'rp cut'). Errorbars, only plotted on two 
of the curves for visibility, assume Poisson statistics. See text for 
discussion. 

by AV and AVs) does not have a strong influence on the dis- 
tribution for either data or mock, except that the 16*'' and 
84"^ percentiles become slightly more noisy as the sample 
size is reduced. There is no strong correlation between veloc- 
ity difi^erence of the satellite-primary pair and the satellite's 
projected distance from the primary. We note that Fig. |4] 
has not been corrected for interlopers, an issue we address 
in Section 133] 

In Fig. [S] we consider the velocity distribution of satel- 
lites averaged over projected radii. Comparing real and red- 
shift space samples, we see that the velocity errors and the 
redshift space distortions tend to broaden the velocity dis- 
tribution. Applying the 'rp cut' to the redshift space mock 
slightly alters the shape of the velocity distribution (the 
number of satellites with smaller velocity difference is re- 
duced more than those with larger velocity difference). The 
dependence of velocity dispersion on project radius seen in 
the mock might not occur in the real Universe. If so the 
incompleteness at small rp has no effect and the 2dFGRS 
satellite samples can be interpreted as sampling the full ve- 
locity distribution within a given radius. However if the real 
satellites are like those in the mock we should restrict our 
comparisons to large scales where the 'rp cut' has no influ- 
ence. 

The velocity distribution of the mocks after applying 
the 'rp cut' is close to the one measured from the 2dFGRS. 
Both sets of satellite catalogues have a velocity distribution 
with an extended tail, with a nearly constant amplitude be- 
yond ~ 800 km s~^. This is to be expected due to contam- 
ination from interlopers and needs to be dealt with when 
analysing redshift space distributions. 

At this stage we split the samples by absolute mag- 
nitude, as our theoretical prejudices, leads us to expect 



Figure 6. Velocity distributions of satellites within a projected 
radius of 400 kpc: similar plot as Fig. \5\ (with same lines 
types showing 2dFGRS data and mock satellite sample in red- 
shift space, with velocity errors and 'rp cut' applied), but split by 
absolute magnitude (indicated in each panel). The y-axis is here 
divided by the total number of satellites in each sample, to facil- 
itate comparisons between samples of different sizes. For clarity, 
we omit errorbars. See text for discussion. 



that the satellite velocity distributions will be more ex- 
tended around brighter primaries. This is exactly what 
we see in the various panels of Fig. |S] The two top pan- 
els show the velocity distributions around the faintest pri- 
maries. For both samples, it seems clear that a velocity limit 
of 1200 km s~^ is probably too large. Although not plotted 
the real space samples do not contain any galaxy with veloc- 
ity AV ^ 450 km s^^). Moreover, the faintest mock sample 
contains a 'lump' of galaxies at AV ~ 900 km s~^. This 
indicates that for faint primaries a smaller value for AVs 
should be chosen, as otherwise contamination from inter- 
lopers will be very strong. The two bottom panels show the 
satellite velocity distributions around the brighter primaries. 
They both show the existence of satellite galaxies with large 
relative velocities, especially the brightest sample for which 
~ 20% (~ 5%) of the satellites have a relative velocity larger 
than 500 km s"'^ (900 km s"'^). Therefore in order to mea- 
sure the velocity dispersion of these systems it is essential 
to sample the full width of the velocity distribution. 

We note that, with the exception for the brightest sam- 
ple, the velocity distributions of the mocks are quite similar 
to those extracted from the 2dFGRS for each absolute mag- 
nitude split sample. For the brightest sub-sample the 2dF- 
GRS velocity distribution is wider than that of the mock. 
We find that this is related to the presence bright 'isolated' 
ellipticals which are found in 2dFGRS sample. We address 
the infiuence of the morphological mix in section 15.21 
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Figure 7. Fitting a Gaussian plus constant to the velocity distribution of satellite galaxies around primaries of three different absolute 
magnitude bins (bright to faint from top to bottom, with magnitude range indicated in the middle column) taken from respectively the 
mock satellite samples in real and redshift space and from the 2dFGRS data (from left to right). The binned satellite velocity distribution 
is shown by the squares with errorbars (assuming Poisson statistic) and the smooth curve the best fitting Gaussian plus constant, as 
determined by solving Eq. |4] The binning used is regular and of width 3/4 CTbeati which, given in each panel, is the best velocity dispersion 
estimate, expressed in km s~^. Intervals without bins represent velocity intervals without any satellites in. The best fit constant, c, is also 
given in s km~^. The inset in each panel show, in dotted, the 1 and 2-a contour levels, in the Sa — logio c plane, of the two parameter fit 
and, in solid, the corresponding one parameter confidence contours. We point out that the velocity dispersion error contours are relative 
to fTbest 8'iid that insets in different panels have different scales. See text for further details. 



4.2 Interlopers and Background Subtraction 

In the previous section, we came across one important con- 
cern for the satellite sample: its contamination by interlop- 
ers. The core of the velocity distributions shown in Figs. [5] 
and [6] are rather well described by Gaussian distributions, 



with some extended tails. These tails are due to interlopers, 
i.e. galaxies which in redshift space just happen to come 
within the selection region, but which in real space are more 
distant and belong to another halo/system. 

Our stacked systems are probing a cylindrical volume 
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in redshift space. Hence to first order we would expect the 
interlopers to be randomly distributed within this cylinder. 
This motivates modelling the velocity distribution of each 
stacked system as the sum of a Gaussian and a constant. 
More precisely we set up a maximum likelihood estimator 
based on the following probability function: 
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where the velocity fitting range is between -vat and vat- 
Hence the maximum likelihood estimator is just the prod- 
uct of the probabilities p(wi) associated with each satellite, 
which can be written as: 
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By maximising this likelihood as function of a and c for 
each sample of stacked primaries, we are able to determine a 
typical velocity dispersion for these systems, together with 
the fraction of interlopers. We note that this approach can 
be applied to any subsample of satellite galaxies. Hence, we 
can, for example, test for a radial dependence of the satellite 
velocity dispersion by just considering satellites in different 
projected radial shells. The a we measure in this way will 
be the underlying velocity dispersion of the stacked satel- 
lite system added in quadrature with the rms error of the 
measured satellite-primary velocity difference. This simple 
subtraction would only be invalid if the measured satellite- 
primary velocity differences were correlated with the mea- 
surements errors and there is no evidence for this in the 
2dFGRS. 

In Fig. [T] we show the results of this method for three 
magnitude bins (bright to faint from top to bottom) taken 
from the mock satellite samples in real (left) and redshift 
space (middle) and from the 2dFGRS data (right hand col- 
umn). Note that we have not subtracted in quadrature the 
'known' velocity measurement error, which for the data and 
the simulation is Vcrr ~ 110 km s~^ (see section [5^ . from 
the fitted velocity dispersion. The inset in each panel shows 
the 1 and 2-a confidence region^: as dotted contours for 
two free parameters (i.e. Ax^ ~ 2.30, 6.17) and as solid 
contours for one free parameter (i.e. Ax^ = 1-0, 4.0). As 
expected, the real space samples (i.e. left most column of 
Fig. [7} are all fit by a pure Gaussian, as in real space we do 
not have the problem of interlopers. The interloper fraction 
averaged over the fitted range is given by 



I(a,c) = 



2cuflt 
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and for the samples in redshift space we find a roughly 
constant value. For galaxies in the range -20.3 > Mb, — 
51ogioh > -20.7, and assuming a cylinder depth of 



Assuming Gaussian errors, the likelihood given by Eq. |4]is dis- 
tributed like a distribution with two degrees of freedom (cr, c). 
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Figure 8. Velocity dispersion of 'stacked' primaries of given ab- 
solute magnitude for various mock samples: in real space with 
and without velocity errors included (filled and open triangles 
respectively) and in redshift space with and without velocity er- 
rors included (filled and open squares respectively) . The errorbars 
plotted are the l-cr errors as obtained from the 2 parameter fits 
on samples taken with cylinder depths of SV = 1200 km s~^, 
for samples brighter than M* , and with cylinder depths of SV = 
600 km s~^, for the fainter samples. The solid line is the median 
satellite velocity dispersion, as estimated from the volume limited 
semi-analytic galaxy catalogue. The two dashed lines shows the 
16*'^ and 84*'' percentiles of the velocity dispersion distribution. 
For samples where the 80 km s~^ velocity uncertainty has been 
included, we subtract from the estimated velocity dispersion the 
total velocity uncertainty, fcrr, in quadrature. See text for discus- 



1200 km s"-', we find interloper fraction of 10.7tlA%- in- 
terestingly the interloper fraction in the mock samples and 
the 2dFGRS data are similar for all the three magnitude 
bins presented in Fig. [7] This is another example of how 
well our mock catalogues mimic the real data. 

Finally, we note that the size of the error on the satellite 
velocity dispersion differs between the mock and 2dFGRS 
samples. The typical l-cr uncertainty on a for the 2dFGRS 
data is approximately 30 km s~^, independent of the best 
fitted velocity dispersion, whereas for the mock samples un- 
certainty is closer to 20 km s^^. This is probably related to 
the greater homogeneity of mock catalogues, which, despite 
their high level of sophistication, do not contain as much 
variety as the real 2dFGRS data (see section [52J- 



5 VELOCITY DISPERSION OF SATELLITE 
SYSTEMS 

We can now estimate the velocity dispersion of satellites 
around stacked primaries for different ranges of absolute 
magnitude within a chosen limiting transverse radial separa- 
tion. The limiting radius needs to be large enough so that the 
composite satellite system contains sufficient satellites, but 
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Figure 9. Velocity dispersion of 'stacked' primaries of given ab- 
solute magnitude for various mock samples: same as Fig. [S] but 
with the 'rp cut' applied. In light grey, we plot, as for reference, 
the same lines as shown in Fig. |8] 

small enough so as to only sample dynamically connected re- 
gions. From theoretical considerations, a limiting radius of 
175 kpc is reasonable for halo masses between 5 x 10^^ 
and 10^^ h~^M0, as it is smaller than their typical virial 
radius, but still large enough to sample a fair fraction of the 
virial volume. Unless otherwise specified, our measurements 
are all done within a projected separation of 175 kpc 
from the primary. 

5.1 Mock Satellite Velocity Dispersion 

In Fig. |8l we plot the estimated velocity dispersion from 
the mock satellite samples in real (triangles) and redshift 
(squares) space, with and without velocity errors included 
(filled and empty symbols respectively). All are in pretty 
good agreement with the distribution of velocity disper- 
sions measured directly from the full semi-analytic simu- 
lation cube, whose median and associated 16 and 84 per- 
centiles are shown by the solid and dashed curves respec- 
tively. We note that when the volume limited simulation 
cube is analysed in this way including or excluding the most 
central satellites has a systematic effect on the measured ve- 
locity dispersion. Discarding all satellites within a projected 
radius of ~ 5 kpc of the primary, results in a velocity 
dispersion which is systematically larger, by 5 to 15%. 

Fig. [5] is like Fig. [S] but with the 'rp cut' applied to the 
samples analysed from both the simulation cube and mocks. 
For this reason, the black solid and dashed curves in Fig. [9] 
are slightly different to those in Fig. |S] (reproduced in grey 
in Fig. [9]). The agreement between the various velocity dis- 
persion estimates is not as good as in Fig. [S] Nevertheless, 
the satellite velocity dispersion inferred after applying the 
'fp cut' agrees within the (typically 50% larger) errors with 
those measured before this cut was applied. In other words. 



Figure 10. Velocity dispersion as measured using the isolation 
criterion as function of associated primary absolute magnitude. 
The data sample is shown for all primary types (filled solid cir- 
cles), but also split by primary galaxy morphology: elliptical-SO 
(open squares) and spiral-irregulars (open circles). The solid and 
dashed lines, same lines as shown in Fig. [9] are plotted here for 
comparison purpose only. See text for discussion. 

the velocity dispersion as measured from satellite galaxies is 
not too sensitive to the innermost spatial distribution. This 
is probably not a big surprise, because of the mixed selection 
effects that come into the isolation criterion. Despite the less 
than perfect agreement between the simulation and mock in 
Fig- El it is important to have demonstrated that the velocity 
dispersion inferred from satellite galaxies is in general agree- 
ment with the median halo velocity dispersion measured in 
the simulation. However, the fact that the agreement is not 
perfect should also not be forgotten when interpreting the 
results of dynamical studies of stacked satellite systems. 

We recall that to arrive at the velocity dispersion es- 
timates plotted as the filled triangles in Figs. [9] and |8] we 
have subtracted in quadrature the mean velocity measure- 
ment error «orr ~ 110 km s~^, which comes from adding 
in quadrature the velocity measurement errors of primaries 
and the satellites. This has the desired effect of producing 
estimates that are all in reasonable agreement with each 
other, but also has the consequence that the fractional error 
on the estimated velocity dispersion is increased by a factor 
(j/(ct^ — Ucrr)'^^^ relative to that of the a that comes from the 
maximum likelihood algorithm. Nevertheless, Figs. [8] and [9] 
show that for nearly the full absolute magnitude range cov- 
ered by our data we are able to recover the underlying ve- 
locity dispersion without any strong biases due to selection 
effects. 

5.2 Satellite Velocity Dispersion from 2dFGRS 

For the 2dFGRS satellite sample, we find it useful to split the 
sample based on whether the primary is an isolated spiral or 
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elliptical galaxy. In Fig. 1101 we plot the velocity dispersions 
for the combined 2dFGRS samples (filled circles) as well for 
the subsets of satellites with spiral (and irregular) primaries 
(open circles) and elliptical-SO primaries (open squares). In 
the mocks, if we infer the morphological types of the pri- 
maries from the bulge-to-disk ratios, then we also find a 
small fraction of elliptical primaries. However, the corre- 
spondence between bulge-to-disk ratio and morphology is 
crude and so we prefer not to split the mock sample in this 
way. Fig. [10] shows that the velocity dispersions from the 
mocks are in very good agreement with those from the over- 
all 2dFGRS sample, but we also see in the 2dFGRS that the 
velocity dispersions of satellites around spiral primaries are 
significantly lower than those of satellites around elliptical 
primaries. 

We also note that for the 2dFGRS the typically error on 
the estimated velocity dispersion is twice as large as in the 
corresponding mock sample. This is most certainly related 
to a combination of the following factors: the mock satellite 
samples are typically two times larger than the two 2dF- 
GRS satellite samples (as the mock catalogues used are not 
split by morphology); due to the limited physics included 
the mock samples are likely to be more statistically homo- 
geneous than the real data; the average number of satellites 
per primary is slightly smaller in the data than in the mocks 
and the smaller total number of satellites in the data will 
lead to more scatter in the estimated velocity dispersion. 

Based on the correspondence between the estimated 
satellite velocity dispersions and those of the underlying 
dark matter haloes found for the mocks in Section 15.11 we 
make the claim that the velocity dispersions of the primaries 
as inferred from the 2dFGRS satellite sample should be a 
reliable tracer of the primaries 'true' velocity dispersions. 
Quantitatively, we expect this method to work for 'isolated' 
galaxies brighter than Mbj — Slogioh ~ —19.0, i.e. galaxies 
which are as bright or brighter than the Milky- Way. 



6 MASS ESTIMATE OF ISOLATED SYSTEMS 

In section O we showed that it is possible, over a range of 
absolute magnitudes, to recover from the satellite velocity 
dispersion the underlying halo velocity dispersion. There- 
fore, it is tempting to go one step further and infer the mass 
which is dynamically enclosed in these systems. There are, 
nevertheless, many issues which needs to be dealt with in 
order to obtain a reliable mass estimate. Two major con- 
cerns are, of course, which mass estimator to use and how 
to choose the radius within which to measure the mass. Our 
approach is to calibrate a mass estimator using the mocks, 
as for them we know the mass of the parent dark matter 
halo of each galaxy. 

6.1 Calibration of mass estimator 

The way we calibrate the halo mass-luminosity relation is 
to measure in the simulation the relation between halo mass 
and halo velocity dispersion and so obtain a way to relate 
the measured velocity dispersion to a halo mass. In Fig. 1111 
we plot the satellite velocity dispersion as measured in the 
volume limited semi-analytic catalogue as function of the as- 
sociated halo mass within 175 kpc, M175. The varying 
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Figure 11. Satellite velocity dispersion in the volume limited 
mock catalogue as function of the associated primary halo mass. 
Both the halo mass and the satellite velocity dispersion are mea- 
sured within a radius of 175 kpc. The different symbols corre- 
spond to isolated systems found in the mock survey and labelled 
as function of primary luminosity (see key). The thick grey and 
black solid lines correspond to the median satellite velocity dis- 
persion as measured from the simulation cube for all primaries 
found in the mocks and all central galaxies in the simulation cube 
respectively. The associated dashed line corresponds to 16^^ and 
84*'^ percentiles of the satellite velocity dispersion distributions. 
For comparison purpose the median dark matter velocity disper- 
sion is also plotted as a thin black line. 

symbol type indicates the absolute magnitudes of the iso- 
lated primaries (see the figure legend) . The galaxy formation 
model used preserves very accurately the mass-luminosity 
hierarchy of central galaxieiQ. By this we mean not only 
that on average bright central galaxies reside in massive ha- 
los and faint central galaxies in the less massive ones, but 
also that the scatter in this whole mass-luminosity relation 
is quite small. This is an essential assumption which needs 
to be accurately satisfied both in the model and the genuine 
data for the method of stacked primary systems to work. 

Mocks created from the same dark matter simulation, 
but with the Cole et al. (2000) galaxy formation model 
have a much larger scatter in the relationship between halo 
mass and central galaxy luminosity and can therefore not 
be used in the calibration process. The scatter present in 
those mocks does not allow satellites of primaries of simi- 
lar luminosity to be stacked, as they can belong to systems 
which are intrinsically too different in mass. Since this work 
began, a whole new suite of semi-analytic galaxy formation 
models with AGN feedback have appeared (e.g. Croton et al. 
(2005), Bower et al. (2006)). These models also have quite a 

^ This is only true for halo masses with circular velocities below 
400 km s~^ , as for halos with larger circular velocities, an artificial 
cut-off in the cooling recipe creates much fainter central galaxies 
in the very largest halos. 
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large scatter in the relationship between halo mass and pri- 
mary luminosity, though less than in the Cole et al. (2000) 
semi-analytic model. We simply restate that by stacking sys- 
tems by primary luminosity and then attempting to infer 
the halo mass, we are implicitly assuming that the scatter 
between luminosity and mass is small and so in calibrating 
such a relation we should use a mock catalogue in which this 
is true. 

Using Springel et al. (2001) semi-analytic model of 
galaxy formation in the full simulation cube, we show in 
Fig. [TT] that using all central galaxies or just those whose 
primaries satisfy the isolation criterion make very little dif- 
ference to the relation satellite velocity dispersion and dark 
matter halo mass. In both cases the medians and the 16**" 
and 84"^ percentiles of the satellite velocity dispersion dis- 
tributions match well over one full magnitude in halo mass. 
They are both well parametrised by a power-law relation of 
the form (T175 oc (A/irs)", with a ranging between ~ 0.42 
and ~ 0.56, depending on which percentile of the velocity 
dispersion distribution one attempts to fit. We note that in 
our mock we do not find any isolated systems residing in 
halos outside the range ~ 5 x lO^^fe"^ Mq ~ 10^^ Mq. 
Hence any mass estimate outside this range is based on as- 
suming the good correspondence found in Fig. [11] around 
isolated halos and central galaxies holds over a larger mass 
range. Finally, we have investigated whether the calibrating 
relation in the mocks is the same for both elliptical and spi- 
ral primaries. For our crude bulge-to-disk ratio assignment of 
morphological type we find the relations are virtually iden- 
tical and so we have opted to use the one overall relation 
present in Fig. [TT]in all cases. 
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Figure 12. Simulation calibrated halo mass, within 175 kpc, 
as function of primary absolute magnitude. The open squares cor- 
responds to the mass of the halos of elliptical galaxies, whereas 
open circles to spiral halo masses. The solid dots corresponds to 
the halo mass of an 'average' galaxy. The errorbars plotted do 
only take into account the error on the measured velocity disper- 
sion. Errors due to uncertainties in the mass-velocity dispersion 
calibration are not taken into account (see text for further de- 
tails). The best fit parameters for the mass-luminosity relation 
given by Eq. |6]are given for spirals and ellipticals, and plotted as 
solid lines. 



6.2 Mass estimates for 2dFGRS primaries 

We now use the relation between the median velocity disper- 
sion of all central galaxies and halo mass, shown in Fig. 1111 
to directly convert the measured satellite velocity disper- 
sions into estimates of halo masses. In Fig. 1121 we plot the 
estimated masses for primaries split by morphological type. 
As is to be expected, we find in both cases that the inferred 
halo mass increases steadily with the absolute magnitude of 
the primary. With this mass calibration, we find that ellip- 
tical galaxies live in halos which are typically 3 to 10 times 
more massive than spirals of similar bj brightness. 

The scatter seen in Fig. [TT] indicates that the uncer- 
tainty in the mass calibration is large. For a given halo mass, 
the scatter in velocity dispersion is around ~ 25 to ~ 40%. 
Hence, considerable care should be taken when using this 
relation to infer halo mass from the measured velocity dis- 
persion. In Fig. 1121 we are assuming that the satellite ve- 
locity dispersion measurement inferred from the 2dFGRS is 
in good agreement with the median of the 'true' satellite 
velocity dispersion. In Fig. [S] we showed that for the mocks 
the corresponding agreement is good, but not perfect. For 
the mocks we can use the spread between 16"^ and 84*'^ 
percentiles of the distribution as a guide to the uncertainty 
in this calibration and hence on the the systematic uncer- 
tainty in the calibration procedure. For an individual object, 
a systematic shift of ~ 75% in the estimated mass is entirely 
acceptable, as this shift corresponds to the 1-a dispersion on 
the calibration relation. 



Using this calibration, we find the power-law fit relation 

between the halo mass within 175 h^^ kpc. Mi 75, and the 
6j-band luminosity of the primary, Lbj. This can be written 
equivalently as 

logio [Mi75/(h-'M0)] = a logio [Lb,,/L.] p, (7) 

where L» is the characteristic luminosity as given by the 
2dFGRS Schechter luminosity function estimate of Norberg 
et al. (2002), and (3 = logio[M./(h~^MQ)], with M. the 
dynamical mass of an L« galaxy. For elliptical galaxies, we 
observe a nearly linear relation between halo mass and lu- 
minosity, with a ~ 1.0 ± 0.4 and /? ~ 12.6 ± 0.2, i.e. a 
nearly luminosity independent mass-to-light ratio. This is 
in stark contrast with spiral galaxies, for which the relation 
between halo mass and primary luminosity is much steeper, 
with a ~ 1.6 ± 0.6 and [3 ~ 11.8 ± 0.3. The errors on these 
best fit parameters are, as already mentioned, rather sub- 
stantial. 

We note that the existence of a difference between the 
scaling relations for ellipticals and spirals is independent of 
the calibration used. The calibrating relation between halo 
mass and satellite velocity dispersion was assumed to be in- 
dependent of galaxy morphology and hence the difference in 
the mass-luminosity scaling relations is due entirely to the 
differences in the measured luminosity - velocity dispersion 
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relations. However the magnitude of this difference does de- 
pend on the halo mass-velocity dispersion calibration used. 

Finally, for a spiral galaxy like the Milky- Way we es- 
timate the halo mass within 175 kpc to be approxima- 
tively 3.5I2 '5 ^ 10^^ h~^M0, whereas for an elliptical galaxy 
of similar bj brightness we estimate its halo mass to be 
nearly eight times larger. We note the large statistical error 
associated to this mass estimate, which do not include any 
systematical error on the mass calibration relation used. 



6.3 Comparison with other mass estimates 

This statistical Milky- Way mass estimate needs to be com- 
pared with other Milky- Way mass estimate, using com- 
pletely different techniques. When comparing with tech- 
niques independent of the value of the Hubble constant, 
we assume here Ho = 70 kms~^ Mpc~^, implying our es- 
timate of the Milky-Way mass within 250 kpc becomes 
5.OI3 Q X 10^^ Mq. Assuming a simple scaling relation for 
the mass enclosed at large radii (e.g. M(r) ~ r, isother- 
mal sphere), we estimate the mass within 100 kpc to be 



-1.2 ^ ^'^ "^Q- 
Prom the dynamics of the Magellanic Clouds and the as- 
sociated stream, Lin, Jones, & Klemola (1995) estimate the 
mass of the Milky Way inside 100 kpc to (5.5±1) lO" Mq. 
From the escape velocity and motions of satellite galax- 
ies, Kochanek (1996) estimates the mass of the Galaxy in- 
side 100 kpc to be (5 — 8) x 10^^ Mq. Using more recent 
kinematic information for Galactic satellites and halo ob- 
jects, Sakamoto, Chiba, & Beers (2003) derive an essentially 
model independent Galaxy mass estimate within ~ 50 kpc of 
5.5tl'l X 10" M0, which corresponds to ~ 10^^ Mq within 
100 kpc. This mass, which is nearly twice as large as the 
one found by Lin et al. (1995), seems to be confirmed by 
Bellazzini (2004), who uses the tidal radii of remote globu- 
lar clusters (between 35 and ~ 200 kpc from the Galactic 
Centre) to provide constraints on the mass profile of the 
Milky Way, independently of kinematic data and yielding 
an enclosed mass of ^ 10^^ ^© ~ 90 ^pc. 

All these estimates are at least a factor of two larger 
than our estimate and the more recent Milky- Way mass es- 
timates are as much as a factor 5 larger. There are several 
systematic effects that could be contributing to this differ- 
ence. Firstly, our estimate is a statistical estimate of the 
mean mass for galaxies of a given luminosity. The scatter 
about this mean relation is large and ought to be taken into 
account. For instance, the scatter of the halo mass satellite 
velocity dispersion calibration relation is 25% for a given 
mass, which translated into the mass luminosity relation im- 
plies a typical 75% scatter for a given luminosity. Secondly, 
our mass estimate for the Milky- Way is obtained using a 
power-law fit to the data, which therefore averages data over 
a range of luminosities. If a power-law is not appropriate and 
only the data in a bin centred at the Milky- Way's luminosity 
were used then the statistical error on our estimate would be 
~ 70% larger. Finally, according to the Milky- Way models of 
Kochanek (1996), the assumption of isothermal sphere be- 
tween 100 kpc to 250 kpc may not be appropriate. Following 
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his figure 7, the scaling is closer to M(r) 



implying a 



~ 50% increase in the estimated Milky- Way mass. 

Taking these statistical and systematic issues into ac- 
count, our mass estimate is compatible with those of Lin et 



Figure 13. Comparison of velocity dispersion estimates as func- 
tion of absolute magnitude from recent analyses using similar 
techniques. For Prada et al. (2003), McKay et al. (2002a), Brain- 
erd & Specian (2003) and Brainerd (2005), the quoted errorbars 
are upper limits on the errorbars shown in their respective analy- 
ses. For van den Bosch et al. (2004), we show their final velocity 
dispersions estimates, for their satellite sample optimised so as to 
have, according to mock catalogues constructed using the condi- 
tional luminosity function, the lowest amount of interlopers. Our 
measurements, within 175 kpc, are the same as presented in 
Fig. [To] See section \7\ but also the Appendix [A] [B] and O for 
further comments on this rather busy plot. 



al. (1995) and Kochanek (1996), but still roughly a factor of 
2-3 too small when compared to the recent estimates from 
Sakamoto et al. (2003) and Bellazzini (2004). Clearly our 
statistical method is not optimal for inferring the mass of 
an individual object. A fairer comparison would be to com- 
pare our results with other statistical mass estimators, like 
McKay et al. (2002) did when comparing with weak lensing 
estimates from SDSS (McKay et al. 2001). 



7 COMPARISON TO SIMILAR METHODS 

There have recently been several complementary studies 
that have used similar techniques to stack satellites around 
isolated primaries (McKay et al. 2002; Prada et al. 2003; 
Brainerd & Specian 2003; van den Bosch et al. 2005; Brain- 
erd 2005). The results of these studies, appropriately con- 
verted into Mb J — 51ogioh where necessary, are compared to 
ours in Pig. 1131 At first sight, one sees very large variation 
between the findings of the different authors. However one 
has to be very careful because both the selection criteria and 
method of estimating velocities vary considerably. A more 
detailed comparison between the different analysis methods 
is presented in the Appendices. The tests presented there 
show why some estimates differ, but also raise new concerns 
over other estimates which on face value appear to agree 
when in reality they probably should not. 
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From the analysis presented in the Appendix |^ we 
learn first that the agreement in the results from McKay 
et al. (2002) and Prada et al. (2003) is, as already claimed 
by Prada et al. (2003), spurious. The selection criteria are 
different enough so that when applied to mocks the results 
are expected to be significantly different. It is clear from the 
paper of Brainerd & Specian (2003) that they have not sub- 
tracted the pairwise velocity measurement uncertainty from 
their velocity dispersion estimate, while Brainerd (2005) 
does take it into account. Despite this (see Appendix IB)) . 
the results from Brainerd & Specian (2003) and Brainerd 
(2005) still do not agree with our findings also extracted 
from the 2dFGRS. We believe this is due to some source of 
extra contamination of their satellite sample, but we are un- 
able to reproduce their results. Moreover, we find that the 
difference between our measurements and those of Prada 
et al. (2003) is mostly due to the different selection criteria 
used to define the samples. The larger errors on our 2dFGRS 
velocity dispersion estimates cannot be explained solely by 
the larger redshift measurement errors of the 2dFGRS data. 
Using the data of Prada et al. (2003), we find their quoted 
errors to be approximatively 30-40% smaller than those we 
estimate. 

As a final conclusion from this comparative work, we 
have to point out that the large variety of results present 
in the literature refiects mostly the large variation in the 
proposed methods of both identifying the isolated systems 
and also measuring the satellite velocity dispersions. All the 
proposed methods introduce biases and unless one applies 
the same selection criterion to a set of realistic mock galaxy 
catalogues one cannot quantify these systematic effects and 
the interpretation of the results remain questionable. 



rp ~ 175 kpc. Within the errors, the velocity disper- 
sion appears to be constant within rp ~ 175 kpc and 
rp ~ 375 h'^ kpc. In this range, the satellite velocity dis- 
persion does not depend strongly on the luminosity of the 
primary. 

Our mock catalogues allow us to calibrate the velocity 
dispersion-mass relation for galaxies selected according to 
the isolation criterion of our 2dFGRS sample. Fig. [12] then 
indicates that elliptical galaxies reside in haloes which are 
at least 4 times more massive than spiral galaxies of simi- 
lar bj brightness. Galaxy like the Milky- Way typical reside 
in dark matter halos of mass ~ 3.5l2'i 

10^^ h-^M0 within 

175 h-^ kpc. 

A key assumption in our analysis is that isolated galax- 
ies of similar luminosity reside in halos of similar mass. It is 
this that justifies the stacking procedure. Our semi-analytic 
models allow us to test the validity of this assumption. We 
find that in one of two semi-analytic models that we have 
investigated, the Springel et al. (2001) model, there is very 
little scatter in the relation between central galaxy luminos- 
ity and halo mass. A different semi-analytic model, however, 
that by Gole et al. (2000), predicts considerable scatter in 
this relation and this introduces large errors in the halo prop- 
erties inferred from a stacking analysis. Mocks constructed 
from this model return an increasing satellite velocity dis- 
persion as a function of primary luminosity which, however, 
deviates systematically from the velocity dispersion of the 
host dark halos. The reasons behind this difference in the 
galaxy formation models are not investigated in detail here 
but it serves to illustrate that significant theoretical uncer- 
tainties remain in the kind of analysis that we have presented 
here. 



8 CONCLUSIONS 

We have developed, tested and applied a method to probe 
the properties of extended dark matter haloes around bright 
galaxies. We do this by carefully selecting isolated galaxies 
in the 2dFGRS and using their faint satellites as tracers of 
the gravitational potential. By stacking systems of similar 
primary luminosity to improve the signal-to-noise, we esti- 
mate the satellite velocity dispersion. Realistic mock galaxy 
catalogues, created from cosmological N-body simulations 
populated with a semi-analytic galaxy formation scheme, 
enable us to relate the measured velocity dispersion of the 
satellite system to the velocity dispersion and mass of the 
underlying dark halo. 

Fig. [To] shows evidence for the existence of dark mat- 
ter halos around typical galaxies. Our sample of satellites 
probes the potential well of the primaries out to several hun- 
dred kiloparsecs and demonstrates that the dark halos ex- 
tend many times beyond the optical radius of the primary. 
This is in agreement with the current theoretical picture of 
galaxy formation in a cold dark matter universe (e.g. White 
& Frenk 1991, Kauffmann et al. 1993, Gole et al. 2000). 
The satellite velocity dispersion increases with the luminos- 
ity of the primary and is much larger for elliptical galaxies 
than for spiral galaxies of similar bj luminosity. 

The total extent of the dark halo is not constrained 
by our data. Although, the satellite distribution extends to 
rp ~ 375 kpc, most of the signal comes from within 
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APPENDIX A: COMPARISON WITH 
SATELLITE VELOCITY DISPERSION 
MEASUREMENTS FROM SDSS DATA 

The studies of McKay et al. (2002) and Prada et al. (2003) 
are based on different releases of SDSS data. In Fig. 1131 
their results are shown by dot-connected open squares and 
by long-dashed connected filled squares respectively. At first 
sight their results appear in rather good agreement with 
each other, and rather different to ours, which are shown by 
large filled circles (all primaries), large open circles (spiral 
primaries) and large open squares (elliptical/SO primaries). 
However, it is essential to consider the differences expected 
due to the differing satellite selection criteria. 

Al McKay et al.: SrUs = 1.5 &z Sn = 2 h'^Mpc 

McKay et al. (2002) (and also Prada et al. (2003) with 
their sample ^Q) use a much less stringent isolation criterion 
on the neighbourhood of the primary. Our requirement for 
the primary to be at least 8 times more luminous than any 
of its satellites is relaxed to just 4 times brighter. We recall 
that the role of this constraint is to avoid including satellite 

^ This sample is presented in their paper, but is not the sample 
for which they derive their main results. 
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Figure Al. Comparison of velocity dispersion estimates as func- 
tion of absolute magnitude, adopting the same criterion as McKay 
et al. (2002) and Prada et al. (2003) in their sample 3. The la- 
belling is the same as in Fig. 1131 with the exception that we 
have added our mock catalogue results (open triangles and filled 
squares for real and redshift space respectively) for this new iso- 
lation criterion, where the velocity dispersion is measured for 
satellites between 15 and 250 kpc using a cylinder depth of 
500 km s~^. The filled large points correspond to our measure- 
ment done on sample 3 of Prada et al. (2003). Note that the 
connected filled squares from Prada et al. (2003) are estimates 
using the selection criterion discussed in Fig. IA2I 

systems in which the potential well and hence the dynamics 
of the satellites is not dominated by the the primary. To 
avoid this being the case it is necessary to ensure that all 
the satellites are much less massive than the primary. The 
recent mass-to-light ratio measurements of Eke et al. (2004) 
indicate that for around 2 L* , which is the luminosity of our 
brighter primaries, a factor of 8 in luminosity corresponds 
to factor 10 in mass. However for lower luminosity primaries 
the corresponding mass factor is smaller and this argues for 
having a large luminosity difference between satellite and 
primary. 

Another difference between our standard isolation cri- 
terion and the one used by McKay et al. is the much larger 
outer exclusion radius. Our requirement is to have all galax- 
ies within 1 h~^Mpc be at least 0.8 magnitudes fainter than 
the primary, whereas McKay et al. require for the same mag- 
nitude difference a distance of 2 h~^Mpc, which is certainly 
a more restrictive criterion. 

Therefore in order to make the appropriate comparison, 
we need to apply the same selection criterion as McKay et 
al. (2002). Fig. I All compares the results of McKay et al. 
(2002) with new estimates we have made from the 2dFGRS 
data after adopting their selection criterion. Also shown in 
Fig. lAll are velocity dispersion estimates we have made us- 
ing sample 3 of Prada et al. (2003). Using the same selection 
criterion, we now find results which are in better agreement 
with the SDSS estimates, but with larger errorbars. There 
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are two reasons why are our errors are larger. Firstly, the 
much larger velocity uncertainty of 2dFGRS redshift mea- 
surements (typically 85 km s~^ on an individual galaxy for 
the 2dFGRS, compared to less than 30 km s~^ for SDSS) 
causes the errors from the maximum likelihood method to be 
increased by 50 %, for systems with a 'true' satellite velocity 
dispersion of 100 km s~^, instead of just 7 % in the case of 
SDSS. Secondly, the errors quoted by Prada et al. (2002) 
are intrinsically smaller than those we obtain when perform- 
ing our maximum likelihood estimation on their sample 2. 
This result can be partially explained if Prada et al. have 
assumed Poisson statistics to derive their quoted errorbars. 

Comparing the results of applying this relaxed selection 
criterion to the mock catalogues, with the distribution of the 
underlying satellite velocity dispersions shown by the heavy 
solid and dashed lines in Fig. lAll we see that there is a large 
bias. The velocity dispersion recovered from the mocks lie 



outside the 16*'' to 



percentile band of the 'true' un- 



derlying galaxy velocity distribution. In our analysis of the 
samples with relaxed isolation criterion we also note that for 
both the 2dFGRS data and the mocks the stacked satellite 
velocity distributions are no longer well fit by the 'Gaus- 
sian plus a constant' model. The samples look much more 
like the real space velocity distributions shown in Fig. [3 i.e. 
there is just an upper limit for the constant. In other words, 
with these satellite samples, fitting the 2dFGRS data with a 
'Gaussian plus a constant' is not appropriate. On the other 
hand, the satellite velocity distribution of sample 3 of Prada 
et al. is well fit by a 'Gaussian plus constant'. This indicates 
some clear difference between the two satellite samples. 

A2 Prada et al.: Smi, = 5ms = 2.0 & 
Srt = 0.5 h-^Mpc 

The selection criterion proposed by Prada et al. (2003) are 
more conservative than those of McKay et al. (2002), as 
they require a factor of ~ 6 in luminosity between the pri- 
mary and any of the surrounding galaxies. However, for the 
spatial isolation of the primary, Prada et al. are much less 
restrictive with an outer isolation radius, Sri,, of 0.5 ft^^Mpc, 
compared to 2 /i~^Mpc for McKay et al. (2002). Compared 
to our choice, this criterion is less restrictive, especially as 
we require within Srg = 0.5 /i~^Mpc a factor of 8 in lumi- 
nosity difference, and within Srt = 1 /i~^Mpc the primary 
to be at least twice as luminous as any other galaxy. How- 
ever, this is not the only difference between the two SDSS 
satellite analyses. The biggest difference is that Prada et al. 
fit for the satellite velocity dispersion and interloper fraction 
separately in each bin of projected radius. Thus they allow 
both the interloper fraction and velocity dispersion to vary 
with projected radius. 

Including this extra freedom in their velocity disper- 
sion calculations gives rise, as we show in Fig. IA2I to a 
substantial difference in the velocity dispersion estimates, 
compared to what is found in Fig. lAll using the analysis 
method proposed by McKay et al. (2002). From Fig. IX2l 
we conclude that satellite velocity dispersion estimates from 
2dFGRS and SDSS (using sample 2 of Prada et al. 2003) 
agree with each other and that the criterion and method 
used seem to be able to recover (within its large uncertain- 
ties) the underlying velocity dispersion, which is not the case 
for the method of McKay et al. (2002). Finally we are able. 
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Figure A2. Comparison of velocity dispersion estimates as func- 
tion of absolute magnitude, adopting the same criterion as used 
by Prada et al. (2003). The labelling is the same as in Fig. 1131 
with the exception that we have added our mock catalogue re- 
sults (open triangles and filled squares for real and redshift space 
respectively) for this new isolation criterion, where the velocity 
dispersion is measured for satellite between 15 and 90 kpc 
using a cylinder depth is 500 km s~^. The filled large points cor- 
respond to our measurement done on sample 2 of Prada et al. 
(2003). Note that the dot-connected open squares from McKay 
et al. (2002) are obtained with the selection criterion discussed 
in Fig. 



but with larger errors, to recover the velocity dispersion es- 
timates found by Prada et al. (2003). Only for the brightest 
galaxies is there any indication of a discrepancy, and this is 
mainly related to error estimates from Prada et al. (2003) 
which for those samples are clearly much smaller than we 
find. 

Regarding the statistics of the isolated satellite systems, 
we note that adopting the selection criterion of Prada et al. 
(2003) causes the number of systems with at least 6 satellites 
within 250 h^^ kpc to increase drastically. This statistic is 
similar to the one we used before, i.e. at least 9 satellites 
within 400 kpc, to identify and remove small groups 
from our sample. Table lAll compares these statistics for the 
2dFGRS data for our original selection criterion, those of 
Prada et al. (2003) and those of McKay et al. (2002). Typ- 
ically the Prada et al. (2003) selection criterion results in 
twice as many 'large systems' as for our standard selection 
criterion, independent of the value of Nvioi and AVs used. 
Hence the sample of primaries might not be as well isolated 
as one hopes. Interestingly, as long as Nvioi = 0, the selection 
criterion of McKay et al. are very comparable to ours as re- 
gards the number of 'large systems' identified. Nevertheless, 
we note that in all cases the fraction of satellite galaxies that 
are in these 'large systems' is always small. 
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Table Al. Properties of the combined NGP & SGP satellite samples around bright galaxies, for two values of iVvioi and for different 
selection criteria, as indicated by Srus and AVs. The McKay et al. selection criterion correspond to Srus = 1.5, Prada et al. to &ms = 2.0 
and ours to Sms = 2.2. All satellites are within 250 kpc from the primary. The columns labelled 'large systems' indicate how many 
primaries and satellites are found in systems with 6 or more satellites within 250 kpc. The last 3 columns show the total number 
of satellites within two cylindrical shells, where Vp is expressed in kpc. These numbers include the satellites belonging to 'large 
systems'. 
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A3 Results from mocks 

With the comparisons of the two previous subsections, we 
have not been able to really quantify the difference between 
the proposed isolation criteria, even though our results hint 
towards the fact that the isolation criterion proposed by 
McKay et al. is probably the least appropriate of the three 
considered and that the one by Prada et al. seems to identify 
slightly more 'large systems' than ours. Therefore, in this 
section, we use the mocks and address the issue of the radial 
dependence of the satellite velocity dispersion, which Prada 
et al. found clear evidence for in their data. 

In Fig. IA3l we show a comparison between the different 
selection criteria used by McKay et al. (left panel) , Prada et 
al. (central panel) and us (right panel). The aim is to show 
how well each of these different selection criteria succeed in 
recovering the 'underlying' satellite velocity dispersion, for 
which the median is given by the solid lines, and the dotted 
lines represent the 16**^ and 84**^ percentiles of the satellite 
velocity distribution as measured from the simulation cube. 
For comparison purposes, we have chosen to present in each 
panel the results within two different projected radii, as in 
Prada et al. (2003): 15 < rp/h'^ kpc < 90 and 175 < 
rp//i-i kpc < 250. 

First of all we note that the different criteria are not 
all as successful in recovering the underlying satellite veloc- 
ity dispersion. Clearly the one proposed by McKay et al. 
(i.e. left panel of Fig. IA3[I is the least successful, as it sys- 
tematically gives an underestimate of the satellite velocity 
dispersion. This is especially true for the outer radial bin, 
for which the measured satellite velocity dispersion is barely 
within the 16"' percentile of the underlying satellite velocity 
dispersion distribution. Moreover we note that for the mocks 
there are not enough faint primary systems for which this 
satellite criterion is satisfied, explaining why no velocity dis- 
persion measurement is given for primaries fainter or equal 
to M* . Interestingly, this is not true for the real data, for 
which there are still several faint primary satellite systems 
satisfying the isolation criterion. With larger number statis- 



tics, this could be a potential way of constraining certain 
galaxy formation models, an approach followed up by van 
den Bosch et al. (2005). 

The middle and right panels of Fig. IA3l look more sim- 
ilar, in the sense that they both recover, within the errors, 
the 'underlying' satellite velocity dispersion, for both pro- 
jected radial bins. However, three small differences can be 
noted. Firstly, the velocity dispersion errors around bright 
primaries are much larger using the isolation criterion of 
Prada et al. (2003). We note that the consistency between 
expectation and measurement is only reached due to large 
non-symmetrical errors. Secondly, the median satellite ve- 
locity dispersion in the outer projected radial bin is less ac- 
curately recovered with the isolation criterion of Prada et al. 
(2003). Indeed, both real and redshift space measurements 
have a tendency to predict a much flater luminosity- velocity 
dispersion relation than the underlying one. Thirdly, the real 
space measurements obtained with the isolation criterion of 
Prada et al. (2003) severely underestimate the true underly- 
ing velocity dispersion, especially for the brighter primaries. 

Finally we note that, for our mock catalogue, there is 
virtually no radial dependence of the satellite velocity dis- 
persion on the choice of radial shell (the distributions in- 
dicated by the bold and shaded lines are all in very good 
agreement with each other) . It is worth pointing out that of 
the three proposed isolation criteria, only the one proposed 
by Prada et al. hints at a radial dependence of the satellite 
velocity dispersion. The trend is definitively weak and would 
be insignificant if if the errors happen to be underestimated, 
or even just assumed to be symmetrical. The fact that this 
can happen is a potential worry for the claim made by Prada 
et al. (2003) for a radial dependence of the satellite velocity 
dispersion. 



A4 Conclusions for SDSS comparison 

We now have to address the question of whether our claim, 
that the isolation criteria of Prada et al. and McKay et al. 




Figure A3. Mock results using the isolation criteria of McKay et al. (left panel), of Prada et al. (centre), and of this paper (right 
panel). The symbols with errorbars correspond to the real (open triangle) and redshift (filled square) space satellite velocity dispersion 
measurements around mock primaries identified using different selection criteria. The lines correspond to the median (solid) and IG*'' 
and 84^^^ percentiles (dashed) of the underlying primary satellite velocity dispersion. The different shadings correspond to two different 
cylindrical shells, whose ranges are given in the panel. 



are more relaxed than our standard ones, is consistent with 
the differences found between the estimated velocity dis- 
persions. Looking at the systematic difference between ve- 
locity dispersion measurements made with the two criteria, 
what we find is quite counter intuitive. One would proba- 
bly expect, that relaxing the isolation criterion would result 
in measuring larger satellite velocity dispersions, whereas 
the opposite is found. This can probably related to a third 
difference between our nominal satellite selection criterion 
and the one used by Prada et al. : the length of the cylin- 
der, within which satellite galaxies need to reside in order 
to be considered in the velocity dispersion estimate. Using 
Monte-Carlo realizations of satellite sample^ drawn from 
a 'Gaussian plus a constant' velocity distribution, we see 
that if one uses exactly the same velocity criterion as in 
Prada et al. (i.e. AVs = 500 km s~^), one starts to sys- 
tematically underestimate the velocity dispersion of systems 
with intrinsic velocity dispersion larger than 180 km s~^. At 
the same time, one systematically overestimate the back- 
ground for those systems. Taking into account the velocity 
errors, this translates, in the case of the 2dFGRS, to systems 
which are best fit by a velocity dispersion of around or larger 
than ~ 210 km s""*^, corresponding to all primaries slightly 
brighter than M* (see e.g. Fig. [T]). On those grounds, we 
motivate therefore the use of a deeper cylinder than Prada 
et al. in order to accurately measure the velocity dispersion 
of slightly larger systems. Our Monte-Carlo approach shows 
that we could have adopted a limiting velocity difference of 
900 km s~^ for spiral primaries, whereas there is a need, for 
elliptical primaries, to go out to ~ 1200 km s~^, in order to 
appropriately sample the satellite velocity distribution. 



° i.e. with same number of systems as in the 2dFGRS, with sim- 
ilar distribution of satellites per system, and with the inclusion 
of an intrinsic maximal 10 to 20 % variation in the background 
and in the underlying velocity dispersion, so as to mimic to some 
extent that all systems are not exactly identical. 



Therefore, most of the discrepancy between our anal- 
ysis and the two using SDSS data reside in the different 
selection criteria used. From our results using the mock cata- 
logues, from our Monte-Carlo simulations and from statistics 
of large systems discovered by the isolation criteria, there is 
a hint in the direction that the relaxed isolation criteria used 
by McKay et al. and Prada et al. are not as appropriate for 
finding dynamically isolated systems as our more stringent 
isolation criterion. Nevertheless we have to point out that 
the data is not yet good enough to be able to fully discrim- 
inate between the methods chosen. 

Finally, a closer inspection of the two SDSS works shows 
that their findings are slightly different, something which 
was already pointed out in the analysis of Prada et al. 
(2003). Indeed the outer radii within which the velocity dis- 
persions are measured are very different. As Prada et al. 
claim a strong dependence of the satellite velocity disper- 
sion on radius, the agreement seen in Fig. [T3]is not as good 
as it looks. On this last point, we would like to add, that we 
are not able, with our isolation criterion applied to the 2dF- 
GRS, to detect such a signal. We know that the behaviour 
for ellipticals and spirals is rather different as function of lu- 
minosity, and therefore it could be legitimate to ask whether 
the effect seen by Prada et al. as function of projected radius 
could be due to a change in their sample mix as function of 
luminosity. Indeed with our findings, for galaxies of similar 
brightness, ellipticals will reside in much larger haloes than 
spirals. Hence stacking galaxies together irrespective of their 
morphological type, as done by Prada et al. (2003), could 
give rise to a velocity dispersion which depends on the ra- 
dius within which it is measured. With our samples we are 
not able to reliably examine this issue, as it requires the 
samples to be split by morphological type, luminosity and 
projected radius. The only conclusion we can draw from our 
samples is that we observe a trend indicating that satellite 
velocity dispersion measurements of galaxies residing in the 
range 175 < Vp / kpc < 375 do not contain much infor- 
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mation. This is in perfect agreement with the fact that the 
meeisured satelhte velocity dispersions within 375 kpc 
are identical, to within the errors, to those measured within 
175 ft-i kpc. 



APPENDIX B: COMPARISON WITH 
BRAINERD & SPECIAN AND BRAINERD 

Regarding a comparison with the measurements of Brainerd 
& Specian (2003), for which in Fig. [13] ellipticals and spirals 
are shown by dashed-connected filled and open pentagons 
respectively, we first need to point out that they have used 
a similar selection criterion to the one proposed in Prada et 
al. (2003). For that reason, we expect, as explained above, 
to find differences between their results and our standard 
ones. However, like we did for Prada et al., we should be 
able to recover their results by assuming the same selection 
criterion. 

For their sample of satellites around spiral primaries, it 
is impossible to recover their results for the following rea- 
sons. Firstly, they have forgotten to subtract the rms veloc- 
ity measurement errors in quadrature, which in their case are 
of the same order of magnitude as ours, i.e. ~ 110 km s~^ 
(as they use data from the 2dFGRS 100k release) . Secondly, 
in some way, the isolation criterion they have applied have 
to be wrong, as it is impossible to understand how they ini- 
tially find an isolated system with more than 605 satellites. 
This is even larger than the largest galaxy cluster found in 
the complete 2dFGRS by Eke et al. (2004). In all likelihood 
they must have forgotten to deal with effects due to 2dF- 
GRS 100k window function, which is extremely patchy, and 
hence very unsuitable for this type of study. Similarly, they 
have probably not used the full photometric input catalogue 
to reject systems which were not fully observed. For these 
reasons, their velocity dispersion measurements around spi- 
ral galaxies, from data which is a subsample of what we used 
in this analysis, are strongly erroneous. 

Regarding their subsamples of elliptical primaries, there 
is no reason not to believe that they are affected by the same 
problems as their sample of spiral galaxies. Ifowever, due to 
the fact that the intrinsic velocity dispersion of those sys- 
tems is much larger, forgetting to subtract in quadrature the 
velocity errors does not infiuence the results by more than 
10 to 15% systematically. On the other hand, we believe, 
due to their problems with the isolation criterion, that the 
errors they quote on the satellite velocity dispersion around 
elliptical primaries is probably underestimated, and that the 
very strong trend with luminosity is too large. 

Finally, we note that we have not, at all, been able to 
reproduce their results using either the full 2dFGRS sample 
nor the 2dFGRS 100k release sample. We suspect therefore 
some of the above mentioned problems to be the cause of 
these difference. 

Recently, Brainerd (2005) made a new satellite analysis 
using the full 2dFGRS survey. Like for the Brainerd & Spe- 
cian (2003) work, we are unable to reproduce in detail their 
findings, especially for faint primaries. 



APPENDIX C: COMPARISON WITH VAN DEN 
BOSCH ET AL. 

Finally, the results from van den Bosch et al. (2004), shown 
by dot-connected filled triangles in Fig. 1131 clearly show 
that, with an isolation criterion that is too relaxed, the pro- 
posed method no longer finds a majority of dynamically iso- 
lated systems. Hence their criteria, as already discussed in 
their paper, were not intended and should not used for se- 
lecting systems for dynamical studies. 

We note that the analytic method presented in van den 
Bosch et al. (2004) for the luminosity - velocity dispersion 
relation, derived for their conditional luminosity function 
model, is in very good agreement with our measurements 
of that same relation from our semi-analytic mocks. More- 
over, applying our isolation criteria to van den Bosch et 
al. (2004)s conditional luminosity function mocks, we re- 
cover the underlying luminosity - velocity dispersion rela- 
tion to great accuracy. This is a very strong consistency test 
for two completely different sets of models, constructed and 
constrained by different mechanisms. 



